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We develop a first quantization description of fractional Chern insulators that is the dual of the 
conventional fractional quantum Hall (FQH) problem, with the roles of position and momentum 
interchanged. In this picture, FQH states are described by anisotropic FQH liquids forming in 
momentum-space Landau levels in a fluctuating magnetic held. The fundamental quantum geometry 
of the problem emerges from the interplay of single-body and interaction metrics, both of which act 
as momentum-space duals of the geometrical picture of the anisotropic FQH effect. We then present 
a novel broad class of ideal Chern insulator lattice models that act as duals of the isotropic FQH 
effect. The interacting problem is well-captured by Haldane pseudopotentials and affords a detailed 
microscopic understanding of the interplay of interactions and non-trivial quantum geometry. 
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The effects of topology and quantum geometry in con¬ 
densed matter physics have recently garnered immense 
attention. With topological insulators and the quantum 
anomalous Hall effect constituting the well-understood 
case of non-interacting or weakly-correlated electron sys¬ 
tems [1 3], the recent theoretical discovery of the frac¬ 
tional quantum Hall (FQH) effect in nearly-flat bands 
with non-trivial topology [4-8] poses deep questions re¬ 
garding the confluence of strong interactions and non¬ 
trivial quantum geometry. On one hand, an experimen¬ 
tal realization of such a fractional Chern insulator (FCI) 
could conceivably push relevant energy scales by an order 
of magnitude, paving the way to robust FQH signatures 
[9-14]. On the other hand, the disparity of conventional 
Landau levels and fiat bands with non-zero Chern num¬ 
ber C suggests a rich playground to realize novel states 
with topological order that cannot be attained in a con¬ 
ventional electron gas in a magnetic held [15, 16], while 
simultaneously presenting a profound challenge to under¬ 
stand the underlying microscopies of strong interactions. 

Most of our current understanding regarding FCIs 
stems from exact diagonalization (ED) of small clusters 
for C = 1 [17-30] and C > 1 [31-35], mutatis-mutandis 
mappings of the Hilbert space of flat Chern bands to the 
lowest Landau level (LLL) [15, 36-40] or vice versa [16], 
and approximate long-wavelength projected density alge¬ 
bra [41-48]. The latter approaches however treat exclu¬ 
sively the universal long-wavelength continuum limit of 
the FQH problem, whereas the presence and relevance 
of the lattice manifests itself in the short-wavelength 
physics. This conundrum is highlighted by the zoo of FCI 
lattice models established so far, which display strongly 
varying proclivities to host stable FQH phases that do not 
correlate well with simple measures such as “flatness” of 
band dispersion and Berry curvature. 

At its heart, the theoretical challenge stems from the 
fact that the immense success in describing the micro¬ 
scopies of the conventional FQH effect resists a simple 
description in second quantization [49] that is essential 


to describe interacting electrons on the lattice [50-55]. A 
resolution is crucial to provide a foundation for studies 
of non-Abelian phases [56-58] and to provide microscopic 
insight that can ultimately drive experimental discovery. 

In this work, we develop a first quantization description 
of FQH states in FCIs that leads to an effective Hamil¬ 
tonian that is the dual of the usual FQH problem, with 
the roles of position and momentum interchanged. In 
this picture, FCI analogues of FQH states are described 
by anisotropic FQH liquids forming in momentum-space 
Landau levels in a fluctuating magnetic field. The chal¬ 
lenge of understanding FQH states in FCIs reduces to 
a variational problem of determining the deformation of 
the guiding-center orbitals due to the presence of the lat¬ 
tice, in analogy to Haldane’s geometrical picture of the 
anisotropic FQH effect (FQHE) [59-61]. Guided by these 
insights, we then present and provide examples of a broad 
class of ideal FCI host lattice models that constitute FCI 
analogs of the isotropic FQHE. These models afford a par¬ 
ticularly simple description of the interacting low-energy 
dynamics, acting as FQH parent Hamiltonians with emer¬ 
gent guiding center and SU(C) symmetry. The effects of 
quantum geometry and Berry curvature fluctuations are 
analyzed in terms of Haldane pseudopotentials. 

Consider a 2D band insulator hosting an isolated 
fractionally-filled flat band with non-zero C, generically 
described by an IV-orbital Bloch Hamiltonian hk. In band 
basis, the flat band of interest is spanned by Bloch states 
|itk) with dispersion hk |ttk) = e |itk}- If the band gap is 
larger than intra-band interactions, then the kinetic en¬ 
ergy is effectively quenched while momentum-dependent 
orbital mixing for |itk) gives rise to a non-trivial quantum 
geometry, expressed by a gauge field and a Riemann met¬ 
ric on CP^ -1 for the Bloch band, the Berry curvature 
Q(k) and Fubini-Study metric g^ v (k): 

Q(k) = e^dk^Ajyik) A^(k) = -i (u k \d ki/ u k ) (1) 
9u, v { k) = \ (<9 fc(j ttk|[l - K)(wfc|]|dfc„«k) + (p v) (2) 
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Here, A„(k) is the Berry connection, and C = 
A f BZ d 2 k fl(k). We use lattice constants a 0 = h = 1. 

The first task at hand is to describe the guiding- 
center basis. In the case of an isotropic free elec¬ 
tron gas in a magnetic field, Laughlin constructed a se¬ 
ries of incompressible FQH trial wave functions [62] for 
odd-fraction filling factors, from a single-body basis of 
radially-localized symmetric-gauge LLL wave functions 
(r|m) ~ z m e~^ r \ / 4 , which are uniquely determined by 
demanding that they are eigenstates of the angular mo¬ 
mentum operator. While angular momentum does not 
readily translate to the lattice, a key observation is that 
such states are simultaneous eigenstates of a parabolic 
confinement potential V'(r) = ^r 2 projected to the LLL: 
[Em' \m'){m'\] V(r) \m) = A (m + 1) |to). 

A natural way to adapt this construction to an FCI is 
to consider anisotropic confinement on a lattice, as a tool 
to determine the guiding-center basis: 

V(v) = \\x^ v x v (3) 

Here, is a unimodular Galilean metric which a priori 
serves as a variational degree of freedom, constrained to 
retain the discrete rotational symmetries of the host lat¬ 
tice, and r = miai + 77 ^ 2 , 777 . 1,2 G Z indexes the unit 
cell with lattice vectors ai, 2 - Placing V on a L x L lat¬ 
tice via appropriate long-distance regularization [63], the 
low-energy dynamics follow from projection onto the flat 
band with projector P = |k)(itk|, and taking L —7 00 : 

V = PV(r)P r = ^ri^lX + ±rf v g vli keBZ (4) 

where 11^ are momentum-space analogues to the usual 
canonical momentum operators, with p = x,y: 

= - w k + a m ( k) [ri„ ri„] = -*e^ n(k) (5) 

Isotropic case: Physical insight may be gleaned by 
identifying (4) with an electron in a magnetic field but in 
momentum-space; 77 ^" = 5^ is particularly instructive: 

V = \ [-tV k + A(k)] 2 + ^ tr g(k) k G BZ (6) 

Here, fl(k) identifies with the magnetic held, and 

electrons scattering at small momenta pick up Berry 

phase factors in analogy to the Aharonov-Bohm effect. 

For benign ’magnetic Held’ fluctuations over the BZ 

f d 2 k[k B n(k) — V\ 2 /Abz < 1 with inverse magnetic length 

k B = and BZ area Abz , the guiding-center basis 

is therefore described by the well-known Landau levels 

on the torus penetrated by flux C, but in momentum- 

space. These momentum-space Landau levels (MLLs) are 

indexed by two quantum numbers 777., n with eigenspec- 

trum e mn = , x , c (m + 1), where m is the MLL index 
V Abz 

and 77 indexes a C-fold degeneracy per MLL [65]. Iden¬ 
tification of the usual guiding center coordinates is thus 
reversed: the MLL index 777 plays the role of the FCI guid¬ 
ing center index and can be identified with discrete Cn 


rotational symmetry if present, whereas n = 0, ...,C — 1 
acts as a component index for C > 1 [ 66 ]. In real space, 
MLL wave functions are radially-localized (Fig. 1 (a-f)). 
Importantly, (4) does not enter as physical confinement 
in the infinite system, thus A —7 0 can be taken in the 
thermodynamic limit. However, A 7 ^ 0 enters as a proper 
energy scale when considering finite-size droplets [67, 68 ]. 

Anisotropic case: Generically, the confinement metric 
77 ^" can be expressed in terms of a complex vector u> that 
obeys 77 ^" = + u) v oj^. The isotropic limit becomes 

u> = [M] t /v/ 2. Corresponding guiding center operators 

7T = fly w" , 7r f = ^ ri M (7) 

obey commutation relations [if, 7 r^] = H when dk = 0 , 
which fixes a phase freedom u> —» u;e l(p . Substituting 
(7) in (4) yields the confinement Hamiltonian in guiding- 
center language that determines the single-body basis: 

V = A 7T7f t + §A ( 7 - fi) ( 8 ) 

Comparison with Haldane’s construction [59] reveals that 
? 7 ^i/ - g p rec i S ely the FCI momentum-dual of the guiding- 
center metric of the anisotropic FQHE. The root cause 
for this duality can be inferred by noting that the con¬ 
ventional FQHE can in fact be formulated both in po¬ 
sition and momentum representation, with single-body 
dynamics and LLL wave functions form-invariant under 
interchange of complex coordinates x + iy and momenta 
k x + iky. This situation is drastically different in FCIs, 
where the discreteness of the lattice necessitates switch¬ 
ing to momentum space in order to retain a first-quantized 
description in terms of continuous coordinates. 

Conceptually, the challenge of devising a microscopic 
description of FQH states in FCIs reduces to a variational 
problem of determining the deformation of the guiding- 
center orbitals upon placing a FQH liquid on the lattice. 
Given an appropriate choice of 77 ^ (k), any FCI can in 
principle be captured by many-body trial ground states, 
constructed from the single-body eigenstates of ( 8 ): for 
instance, given the lowest MLL wave function 'l'o(k) 
and ladder operators al that generate higher MLLs, the 
Laughlin state at v reads dT ~ ~ aJ-) 1 ^ 1 / ’®'o(k). 

A preferred guiding-center metric can be readily iden¬ 
tified by demanding suppression of a residual dispersive 
term rj^g^ — B in ( 8 ) that delocalizes the MLL basis: 

77 (k) = sj det g(k) g _ 1 (k) (9) 

The dispersion vanishes exactly if and only if 

2 i/ detg(k) = | 0 (k)| ( 10 ) 

This is the condition for an ideal FCI droplet , and is sat¬ 
isfied by every two-band model [63] while placing con¬ 
straints on models with three or more bands. 

Interactions: Insight into stabilization of trial states 
necessitates recasting the many-body problem into first 



quantization. Analogous to the FQHE, density interac¬ 
tions may be generalized via an interaction metric fj^ v \ 

Hi = 1 q PqP-q, Pq = CkaEJ^ ] C kQ (11) 

q o;k 

where 7^ [i^"] = e q ^ dk » enters as a metric-dependent 
momentum-space translation operator. In contrast to 
the FQHE on the plane, translation symmetry and pe¬ 
riodicity in q constrain fj^ to SL(2,Z). Importantly, 
the guiding center metric is still a variational degree 
of freedom if the deviation of (9) and fj^ v is signifi¬ 
cant. To proceed, note that if the ideal droplet con¬ 
dition (10) is satisfied, then an exact operator identity 
[1 — | Uk) <u k |] ui^dk |u k ){«k| = 0 [43] entails that any op¬ 
erator of the form O = A~ (—id^u}^)A + with 

analytic functions A ± can be projected to the flat band 
as POP T = A - (7r)A + (7T^). Consider thus a decomposi¬ 
tion of fj^ + x IJ 'Q' y with complex vector x : if 

X is momentum-independent, then the translation opera¬ 
tor is precisely of this form: Tq = e q+u>>ldk ^e q ~ ul ^ dk ^ with 
q± = QnX 11 — this case is considered in detail be¬ 
low [69]. Note that in the isotropic limit fj^ v ,r) lxl ' = 5^, 
Tq reduces to the conventional translation operator Tq = 
e »q-(-*v k ) w jth q ± = {^q x -j- iq y j/^/2 : while Hj is just the 
usual density interaction. While emphasis has thus far 
been placed on narrowing down to a suitable class of in¬ 
teractions, substantial progress has been made: projected 
to the flat band, the joint dynamics of (3), (11) can now 
be succinctly expressed in guiding-center language: 

H = A A + ^ Vq e ^+(A-*-) e *9-(-!-*]) (12) 

i i<j q 

This Hamiltonian is the central result of this paper - 
it provides a first-quantized description of the low-energy 
dynamics of an ideal FCI in terms of the quantum ge¬ 
ometry r] fJ ’ v ,fj^ 1 ' of the lattice. Its interaction describes 
two-body momentum-space magnetic translations, and 
acts solely on the relative guiding center indices. A key 
consequence is the approximate conservation of center-of- 
mass guiding center, quantified by Berry curvature fluc¬ 
tuations averaged over the BZ with Asz || [ 7 ^, 7 r| m ] || 2 = 
f BZ d 2 k [fi(k) - ^] 2 , where 7T cm/rel = (ft, ± tt 2 )/\/2 
span the two-body problem in the Chern band [70]. Fur¬ 
thermore (12) does not act on the intra-MLL index n, 
stipulating an emergent SU(C) symmetry for C > 1 [71]. 

The physics of the above Hamiltonian can be studied 
via a pseudopotential decomposition of the two-body in¬ 
teraction matrix elements = (mM\ H \m'M') |a=o 

with ?n, M relative and center-of-mass guiding center 
indices, and intra-MLL indices omitted. Approximate 
center-of-mass conservation in ( 12 ) entails that the two- 
body repulsion depends only on the relative coordinate, 
VmmJ ~ Vmm'fiMM'■ Since the guiding center index iden¬ 
tifies with discrete rotational symmetries, it is tempting to 
speak of an emergent continuous rotational symmetry in 
the flat band however, it persists even for anisotropic; 
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FIG. 1. (color online), (a-f) Real-space lattice guiding-center 
wave functions for the ideal isotropic C = 3 model, evalu¬ 
ated from (8). Insets depict analogous conventional LLL wave 
functions, (g) Berry curvature H k (inset) and associated MLL 
spectra for the FCI models discussed in the main text, m 
is the guiding center index; dotted lines indicate ideal flat-fl 
spectrum e m = . The C-fold degeneracy of MLLs reflects 

the C-component basis for C > 1 models. The guiding-center 
structure remains robust in the presence of fluctuations of fl k . 


Cn symmetry of Vq instead constrains relative guiding 
center transitions V m ^ m i to (to — m'jmodN = 0. The 
dominant matrix elements are thus well-captured by Hal¬ 
dane pseudopotentials [72] V m = -fa Em=o Itw-kxu 
which indicate stabilization of FQH trial states. 

The propensity of the MLL basis to lead to a well- 
defined pseudopotential expansion for such models is a 
key advantage of the first-quantized formalism. Treating 
S2(k) fluctuations as a perturbation with ladder operators 
7Tj ->• a,, | m,M) = (d 1 t e i) m (“lm) M |0,0) /Vro!M! leads to 


1mm' — / dq 


Vq(iq-) m ~ m rFl 


m+1 

m-m -\-1 7 


2 q+q- 


k ^ m, y / 2 m ' m m'\/m\{m — m')\ 


(13) 


Here, iF\(-) is the Rummer confluent hypergeometric 
function. The well-known pseudopotentials of the con¬ 
ventional FQHE can be readily recovered for to = m/, 

q± = {q x ± iq v )/V2 With V m = f d 2 q V q L m (^)e~ q2/n . 

Isotropic ideal FCI models: While the focus so far has 
been placed on developing an accurate language for the 
generic anisotropic case, a key follow-up question con¬ 
cerns instead applying above results to find an FCI ana¬ 
log of the isotropic Landau level, particularly favorable for 
FQH phases. Such models indeed exist: the isotropic case 
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FIG. 2 . (color online), (a)-(c) Haldane pseudopotentials for 
the C = 1 three-orbital and C > 1 ideal droplet models 
(on-site and NN interactions). Only odd (even) pseudopo¬ 
tentials determine interactions between same-species fermions 
(bosons). Error bars quantify the residual center-of-mass M 
variation, strongly suppressed for the ideal C — 2 , 3 models 
due to emergent conservation of guiding center, (d) Log plot 
of guiding-center and SU(C ) symmetry breaking in the inter¬ 
acting problem with matrix elements quantified as 

the ratio of the norm of AM = mi + m2 — m3 — tua (AN = 
ui +712 — ri3 — 714 mod C component) symmetry-breaking and 
-preserving two-body matrix elements. The guiding-center in¬ 
dex identifies with C4 symmetry, (e) Role of broken rotational 
symmetry on the lattice, depicted via leading-order V^, m+4 of 
generalized pseudopotentials V^, m * = (m, M\H\m', M). 

r/' 1 ', fj IJV = 5^ v with corresponding Fubini-Study metric 
g xx — 9 yy ; 9 xy = 0, trg(k) = Sl(k) is uniquely satisfied by 
any Bloch state that can be written without normaliza¬ 
tion as a meromorphic function |uk) = | Uk^+iky)- The 
number of poles in the BZ defines C [73]; periodic bound¬ 
ary conditions in k restricts ^Uk^+ik y ) to elliptic functions, 
constrained to C > 2. Skew-anisotropic guiding center 
metrics ensue from distortions of the lattice: for instance, 
a ’squeezed’ FQH liquid with rj xx = l/(,r]yy = ( follows 
from BZ strain deformations with k x —> fk x . 

To illustrate our construction, we consider two multi¬ 
orbital toy models on the square lattice for C = 2, 3 [76], 
with the guiding-center basis of MLLs depicted in fig. 1: 

| s k =2 ) = [ 1, « p{h x + iky) ] T (14) 

j«£ =3 ) =[!>/? p(kx + iky), 7 p'(k x + iky) ] (15) 

Here, p(z) is the Weierstrass elliptic function with periods 
27 t, 27t* and a second-order pole at the T point, and a = 
5.77, ft = —7.64, 7 = 6.73 are band structure parameters, 
chosen to minimize Berry curvature fluctuations. 

The corresponding Bloch Hamiltonian is not unique; 
a possible definition with flat bands is hc(k) = 
1 — /(ufejufc)- In general, hc(k) exhibits 

long-ranged but exponentially-decaying hopping terms 
[73] and acts as an artificial toy lattice model for FCI 
states with topological order, much like the Hubbard 


model truncated to nearest-neighbor (NN) hoppings 
acts as a toy model for strongly-correlated states 
with conventional order. While a classification of 
such ideal FCI host lattice models remains an im¬ 
portant open task, note that simple physical models 
can emerge after truncation of irrelevant hoppings. 
For instance, hc =2 is well-described by a canonical 
d-wave lattice model he = 2(k) = d(k) • <? with d(k) = 
[f(cos k x + cos k y ), t(cos k x — cos k y ), t' sin(fc x ) sin(fc y )] T 
with nearest- and next-nearest-neighbor hoppings. 

Analysis: In addition to the ideal isotropic models 
(14,15), we study both a conventional three-orbital C = 1 
model on the square lattice [50] that does not satisfy 
(10), and an “optimized” variant with flat Berry curva¬ 
ture (var(H k ) « 10 -6 ), obtained via adiabatically adding 
symmetry-preserving hoppings up to 5 th neighbor [63]. 

Fig. 2 depicts the pseudopotential decomposition for 
on-site and NN repulsion, using the MLL basis of Fig. 
1. The C = 2,3 ideal lattice models display emergent 
guiding-center and SU(C) symmetries, manifested in a 
vanishing center-of-mass dependence of pseudopotentials 
(figs. 2b, 2c), and suppression of symmetry-breaking two- 
body interaction terms (fig. 2d). As anticipated from 
(13), on-site repulsion results in a non-zero pseudopoten¬ 
tial only for Vq, acting as an optimal (221)-Halperin state 
[74] parent Hamiltonian for hardcore bosons, whilst not 
stabilizing fermionic FQH states [75]. The latter can be 
remedied by tuning NN repulsion to tune Vi, V 3 . The 
controlled expansion in pseudopotentials is a key merit of 
this construction and highlights that, contrary to common 
perception, the confluence of flat Berry curvature and lo¬ 
cal interaction does not stabilize a fermionic FQH liquid 
in an FCI. Conversely, the non-optimal C = 1 models vi¬ 
olate (10) and do not pin the guiding-center metric to g. 
The guiding-center description (12) is incomplete, broad¬ 
ening the effective interaction range with non-vanishing 
and decaying V m with substantial center-of-mass devia¬ 
tion even for a purely local interaction (fig. 2a) that per¬ 
sists even for uniform Berry curvature. We stress that 
this highlights the shortcomings of long-wavelength limit 
arguments [41-47] in predicting the microscopies of the 
FQHE on the lattice. 

In summary, we introduced a first-quantized descrip¬ 
tion of FCIs, with the FQHE emerging in a picture of 
anisotropic momentum-space Landau levels in a fluctuat¬ 
ing magnetic field. We presented a novel class of ideal 
FCI lattice models as duals of the isotropic FQHE and 
demonstrated their optimality via an expansion of local 
interactions into Haldane pseudopotentials which can be 
determined straightforwardly in first quantization. A pri¬ 
mary goal of this work is to establish a deeper microscopic 
understanding of the stabilization of FQH states in flat 
Chern bands - the resulting interplay of topology and ge¬ 
ometry to determine long- and short-wavelength physics 
on the lattice serves as a natural application of the for¬ 
malism of the anisotropic FQHE. The results presented 
set a foundation for microscopic analysis of non-Abelian 
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phases on the lattice and extension to fractional topolog¬ 
ical insulators. 
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Supplementary Material 


Momentum-Space Landau Levels 


As described in the main text, the guiding-center basis in a FCI can be understood in terms of momentum-space 
Landau levels (MLLs). We seek the low-energy momentum-space dynamics induced by a real-space confinement 
potential V(r) = ^x^g^x^. The confinement metric g^ can be expressed as an N d - fold superposition of reciprocal 
lattice vectors rgid^d", where = mibi + TO2b2 and mi,m 2 € Z. A minimum of three basis vectors d, 
is required to reproduce the three independent components of 77. Regularization of the long-distance behavior of the 
confinement potential on an L x L lattice is not unique; a simple choice reads: 

\ t 2 Nd 

Keg = ^ 7^X2 J2 ^ t 1 _ C ° S ( di ' r / L ^ ( 16 ) 


The regularized form of the confinement potential is used for numerical results presented in the main text. Choosing 
a sufficiently large lattice, one can solve for the guiding-center basis as eigenstates of Ho + K eg , where care must be 
taken to ensure that the band gap of Ho is much larger than the depth of the confinement potential. The lattice 
dimensions are taken to be 160x160 sites for numerical results in this work. As is canonical in numerical studies of 
FCIs, we assume that the energy scale of interactions in the flat band is much greater than its band width and work 
with a modified lattice Hamiltonian Ho = —Ho to eliminate residual dispersion ek hr the flat band. This step is not 
necessary for the ideal models of equations (14), (15) in the main text, as they intrinsically exhibit a perfectly flat 
Chern band. The low-energy dynamics follow from projecting Keg onto the reduced Hilbert space built from the set 
of Bloch states |uk) spanning the flat band, and taking the thermodynamic limit: 


A 


H = lim - Y I k ) ( M k | Keg | u v ) (k' | 

L—>-00 Z z — 


kk'GBZ 


XL 


N d 


= lim — Y l k ) ( u k| Y^ r ! i 

L— >00 4 *■ — ' *■ — ' 


\ ^ d r 

= W Y^ / rf2 k l k )( w k| (-id, ■ Vk ) 2 |uk)(k| 

1 , J BZ 


keBZ 


2 |u k )(k| - 


M k+^/\ k + T 


u k-4i)(k- L 


' N d 


= y / ^ 2 k Y^ d i dV i l k ) 


{-id^-idv) + iA^d v + iA^d^ - ((ukl^^lwk)) 


( k l 


J d 2 k jk) [(— id^ + Afj) g^ u (—id,, + A v ) + rf‘ v g^v] (k| 


(17) 


where A M = — i (uk|d/t|«k) is the Berry connection, O = e IJIJ d tl A v is the Berry curvature, and the Fubini-Study metric 
g^ is defined as: 



(^Wkl^Wk) - (c^Wkluk) (Ukl^Uk) + (M t> y) 


(18) 


The last line of equation 17 uses the fact that the metric tensor g^ v is symmetric, hence g^e^ = 0. In terms of 
operators 11^ = —idy, + Ay, one finally arrives at the first-quantized Hamiltonian (4) described in the main text. 

A smooth gauge choice separates A(k) = -f^[—k y , k x ] T + Afl uct (k) into a topological contribution and fluctu¬ 
ations. C / 0 then stipulates that wave functions acquire twisted boundary conditions V’nm(k + TOibi + 777^2) = 
e~' l T^'- mia2 ~ rn2&i]k il) nm (]i) that constrain the intra-MLL index to 0 .. .C — 1, giving rise to the MLL spectrum described 
in the main text. As a side note, residual dispersion in the Chern band can be understood as a momentum-space 
confinement potential, giving rise to a Darwin-Fock eigenbasis [64] on the torus, with the ratio of confinement A and 
bandwidth entering as a variational degree of freedom when building the single-body basis. 


Ideal Droplet Condition — Two-Band Models 


The ideal droplet condition 2y / clet~g(k) = fl(k) is automatically satisfied by any two-band FCI model. To see this, 
consider a generic two-band Bloch Hamiltonian h(k) = e(k) d(k) • <r parameterized by unit vector d(k) and dispersion 



±e(k), where ef are the Pauli matrices in orbital basis. The Berry curvature and Fubini-Study metric read: 


f2(k) 


9^( k ) 


) d < k > * 


<9d(k) <9d(k) 

X 


dk n 


i ( X d(k)' 

4 ^ dk„ y ’ 


d(k) x 


dd(k) \ 

0 k„ J 


(19) 

( 20 ) 


The determinant of the Fubini-Study metric follows from application of quadruple-product identities and satisfies 
2i/det g(k) = fi(k). 


Haldane Pseudopotentials 


This section discusses the derivation of Haldane pseudopotentials for momentum-space Landau levels. The starting 
point is the two-body interaction in first quantization that follows from equation (12) in the main text. The guiding- 
center basis can be found via numerical solution of the single-body problem. Given the lowest MLL state |*o) and a 


set of ladder operators a, a', 




= 1 generating the m-th MLL | in) = 


1*0 ), one can define two-body states 


M 


,M) = 


V / ro!M! \ V2 


V2 


l*n 


( 21 ) 


The Haldane pseudopotentials follow 


VmM = (m,M 




m, M 


( 22 ) 


and can be evaluated numerically, with the residual center-of-mass dependence studied in the main text. To recover 
the limit of flat Berry curvature, note that the guiding center operators can be written in terms of the ladder operators 
'k i = a.Jks with ks = \j2iz/C the magnetic wave vector for a square lattice, which allows for an algebraic evaluation 
of the two-body matrix elements: 


V mM — I d Q Vq 
m! M' J 


d 2 q V„ 


mM 

1 


V 


ml M' 


= J dqV q 
= J d 2 q V q 

= [ d 2 q V q 


q \Jm\m'\ 
1 


(«i - a 2 ) 


0 iq + (a x -a^)/k B p iq-{a\-a\)/k B 


(d\ - o£) 


0 ) Smm' 


y/mbn 'T 


(v /2i q+/ k B) l {V2iq-/kB) 1 ' 


l\l'\ 


(d) m+l (at) 


m +t 


0) Smm' 


1 ^(V2iq + /k B Y(V2iq_/k B r- m '+\ 

yr - n/~ , Ml - {m + iy.SMM' 


\fm\m !T 


l\(m — m' + l)\ 



( s/2iq-/kB) m - m ' i F x (to + 1; m - m' + 1, 

(m — m')\ 


J MM r 


(23) 


where \Fi(a]b;z) is the Kumrner confluent hypergeometric function. Here q+ = x^g^, 9 - = X m 9m> an d 9 = 9-9+ = 
q/iX^X^Qv While the center of mass guiding center dependence is eliminated, lack of continuous rotational symmetry 
on the lattice does not forbid subdominant scattering between states with different relative guiding center indices to, m! 
unless the parent interaction is purely local Vq —» V. Instead, Cn symmetry can provide weaker selection rules with 
V mrrl > 0 only for m — m! mod N = 0. The conventional Haldane pseudopotentials follow from imposing to = to/ , 
9 ± = {qx ± iq y )/y/2, and expressing the hypergeometric function in terms of a Laguerre polynomial: 


V m = J d 2 q Vq L m 


SI ) e -q 2 /fe| 


(24) 


which is the formula quoted in the main text. 
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Elliptic Function Models 


The main text describes a class of ideal FCI models which can be defined from unnormalized Bloch states |Wfc x +jfc H ) 
written as elliptic functions. To see this, consider a generic N-band model with unnornralized Bloch state written in 
terms of complex momenta 

|ilk) = [<i>i(z,z), (j> 2 {z,z),...,(j) N {z,z)] (25) 


where z, z = k x ±ik y . The corresponding Berry curvature and Fubini-Study metric may be evaluated straight-forwardly: 


N 


° = 2 E 

n>rr 

N 

tr 8' = 2 E 


| tpmdz&n findzfiml ^n^z^r) 


W 




n>m 

N 


| (ftmdztfin (frndztfiml + \ ^m^z^n ^n^z^ml 


\*\ 4 


\<t >\ 4 


_ 0 4*ndz*fim) ^n^z^m) “1“ C.C. 

gxx ~ g yy ~ z 2^ UJ4 

n>m 

E {$mdz4*n $ndz*fim) (4 > mdz4 ) n 4*ndz$m) C.C. 

' 

n>m 

i 


W 4 


det g = - [(tr g) 2 - ( g xx - g y V ) 2 - 4g 2 y 


(26) 

(27) 

(28) 

(29) 

(30) 


One can see by inspection that the conditions 12 = tr g = 2-^/det g and g xx — g yy = g xy = 0 are satisfied if and only 
if dz4>n = 0. This constrains the unnormalized Bloch state to meromorphic functions in k x + ik y . The corresponding 
droplet confinement metric is just the Euclidean metric as described in the main text. Generally, one can consider an 
arbitrary static lattice deformation defined by transformation k —> U • k. As such a transformation preserves both 12 
and det g, 12 = 2 ^/detg remains satisfied, whereas the confinement metric changes to deform the shape of the droplet 
accordingly. 


Three-Orbital C = 1 Model — Flat Berry Curvature Optimization 


To study the effects of Berry curvature fluctuations, we introduced in the main text a C = 1 model with flat Berry 
curvature. This model is adiabatically connected to the three-orbital model by Sun et al. [50] by adding longer-ranged 
hopping terms to minimize Berry curvature fluctuations. We start from the original lattice Hamiltonian 


^ E [ "d.k S x ,k r 'Py. k 

k 


—2tdd (cos k x + cos k y ) + S 2 it p d sin k x 

—2it p d sin k x 2t vv cos k x — 2if cos k y 

—2 itpd sin k y —iA 


2 itpd sin k y 


C d,k 

iA 


C Px,k 

2 tpp cos k y — 2tpp cos k x 


. C Py,k . 


Pj/. k J 

(31) 


for a three-orbital model on the square lattice with d x 2 _ y 2 , p Xl p y orbitals per site. This model exhibits a C 4 rotational 
symmetry and mirror symmetries with symmetry operators 


II 

' 1 0 o' 
0 0 1 

' Rtt/2 Um — 

'10 0' 
0 1 0 


0-10 


0 0-1 


It - K 


(32) 


where R n / 2 , I x and K are 7 t/2 rotation, x-rcflection and complex conjugation operators, respectively. To proceed, we 
write down the most general model Hamiltonian with up to 5 th neighbor hoppings between orbitals that is invariant 
under the symmetry operations above. The parameter vector v of hopping amplitudes serve as the variational parameter 
of the numerical minimization of Berry curvature fluctuations using cost function 



Q(k,v) 



(33) 
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The minimization procedure uses the original 3-orbital model parameters as starting values for v. Furthermore, care 
has to be taken to ensure that the band gap remains open. We quote here the resulting model Hamiltonian that is 
used in the main text: 

hdd hdx hdy 
h* h h 

n dx n xx n xy 

h * h* h 

,L dy "'xy ll, yy 

with 




hdd = 0.203 cos(k x — 2 k y ) — 0.0443 cos(2fc x — 2 k y ) + 2.75 cos (k x — k y ) 

+ 0.203 cos(2 k x — k y ) + 2.75 cos (k x + k y ) + 0.203 cos(2 k x + k y ) + 0.203 cos(fc x + 2 k y ) — 0.0443 cos(2 k x + 2 k y ) 


+ 4.13 cos (k x ) — 0.663 cos(2 k x ) + 4.13 cos (k y ) — 0.663 cos(2 k y ) — 0.76 (35) 

hdx = —(0.0372 — 0.1273?) sin (k x — 2k y ) + (0.0391 — 0.0043?) sin {2k x — 2 k y ) — (0.18 — 1.46?) sin (k x — k y ) 

+ (0.0439 — 0.0177?') sin(2fc x — k y ) + (0.18 + 1.46?) sin (k x + k y ) — (0.0439 + 0.0177?) sin(2 k x + k y ) 

+ (0.0372 + 0.1273?) sin(fc a: + 2 k y ) — (0.0391 + 0.0043?) sin(2 k x + 2 k y ) + 1.01?sin(fc x ) — 0.374? sin(2fc x ) 

+ 3.10sin(fc y ) + 0.0351 sin(2 k y ) (36) 

hd y = (0.0439 + 0.0177?) sin(fc x — 2 k y ) + (0.0391 + 0.0043?) sin(2fc x — 2 k y ) — (0.18 + 1.46?) sin(fc x — k y ) 

— (0.0372 + 0.1273?) sin(2fc x — k y ) — (0.18 — 1.46?) sin(fc x + k y ) — (0.0372 — 0.1273?) sin(2 k x + k y ) 

+ (0.0439 — 0.0177?) sin(fc x + 2 k y ) + (0.0391 — 0.0043?) sin(2 k x + 2 k y ) — 3.10sin(fc x ) — 0.0351 sin(2 k x ) 

+ 1.01? sin (k y ) — 0.374? sin(2fcy) (37) 

h xx = —0.0912 cos(fc x — 2 k y ) + 0.0912 cos(2fc x — 1 k y ) + 0.0912 cos(2 k x + k y ) — 0.0912 cos(/c x + 2 k y ) 

+ 0.156 cos (k x ) + 0.0472 cos(2 k x ) — 0.156 cos (k y ) — 0.0472 cos(2 k y ) (38) 

h xy = (0.0808 + 0.0062?) cos (k x — 2 k y ) + (0.0273 + 0.0191?) cos(2 k x — 2 k y ) + (0.627 — 0.847?) cos (k x — 1 k y ) 

+ (0.0808 + 0.0062?) cos(2fc x — 1 k y ) — (0.627 + 0.847?) cos (k x + k y ) 

— (0.0808 — 0.0062?) cos(2 k x + k y ) — (0.0808 — 0.0062?) cos (k x + 2 k y ) — (0.0273 — 0.0191?) cos(2 k x + 2 k y ) 

— 0.921? cos(fc x ) — 0.0454? cos(2/c x ) — 0.921? cos ( k y ) — 0.0454? cos (2k y ) — 1.07? (39) 

h yy = 0.0912 cos(fc x — 2 ky) — 0.0912 cos(2fc x — 1 k y ) — 0.0912 cos(2 k x + k y ) + 0.0912 cos(fc x + 2 k y ) 

— 0.156 cos (k x ) — 0.0472 cos(2 k x ) + 0.156 cos (k y ) + 0.0472 cos(2 k y ) (40) 


which yields a measure of Berry curvature fluctuations 
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d 2 k f2(k, v) 


1 

277 


= 5.37 x 10 


(41) 



